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A literature surrey has been conducted to assess the state-of-the-art of 
time- temperature dependent elastic-viscoplastic constitutive theories which 
are based on the nnified approach. This class of constitutive theories is 
characterized by the use of kinetic equations and internal variables with 
appropriate evolutionary equations for treating all aspects of inelastic de- 
formation including plasticity, creep, and stress relaxation. The review 
identifies aore than ten such unified theories which are shown to satisfy the 
uniqueness and stability criteria imposed by Drucker's postulate and Ponter's 
inequalities. These theories are coopered on the basis of the types of flow 
law, kinetic equation, evolutionary equation of the internal variables, and 
treatment of temperature dependence. The similarities and differences of 
these theories are first outlined in terms of mathematical formulations and 
then illustrated by comparisons of theoretical calculations with experimental 
results which include monotonic stress— strain curves, cyclic hysteresis loops, 
creep and stress relaxation rates, as well as thermomechanical loops. Numer- 
ical methods used for integrating these stiff time-temperature dependent con- 
stitutive equations are also briefly reviewed. 

INTRODUCTION 

Constitutive theories based on the classical concepts of plasticity and 
creep generally decompose the inelastic strain rate into a time— independent 
plastic strain rate and a time— dependent creep rate with independent constitu- 
tive relations describing plastic and creep behavior. While this approach can 
be rationalized on historical grounds and perhaps on computational conve- 
nience, experimental evidence collected on structural alloys at elevated tem- 
perature indicates inherent time'-dependency and creep/plasticity interactions 
[1] . This suggests that inelastic deformation might be primarily controlled 
by a single overall mechanism and should be treated in a unified manner. 

In recent years, a number of formulations of elastic-viscoplastic con- 
stitutive equations have been presented in the engineering literature. Such 
equations are sometimes referred to as "unified" since inelastic deforma- 
tions are represented and treated by a single kinetic equation and a discrete 
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set of internal variables. In this context, creep, stress relaxation, and 
plastic flow are different manifestations of time-dependent inelastic deforma- 
tions under particular loading conditions with consequently different response 
characteristics. 

There are more than ten unified constitutive theories in the literature. 
These constitutive equations have some common properties and some essential 
differences which have been reviewed recently by Walker [2] . Since then, 
there have been more advances in the development of the unified theories. The 
purpose of this survey is to update Walker's previous work by reviewing the 
state-of-the-art and the numerical integration techniques for these unified 
theories. This survey also serves to identify areas for further model devel- 
opments. 

The unified theories which are reviewed in this survey include those of 
Walker [2], Bodner and Partom [3,4], Miller [5], Krieg, Swearengen and Rhode 
[6], Chaboehe [7], Robinson [8], Hart and co-workers [9], Stouffer and Bodner 
[10], Lee and Zaverl [11], Ghosh [12], and Sagawa and Asada's modification of 
Miller's model [13]. 


GENERAL CHARACTERISTICS OF UNIFIED CONSTITUTIVE EQUATIONS 
FOR ELASTIC-VISCOPLASTIC MATERIALS 


Constitutive equations for elastic-viscoplastic material could be formu- 
lated either with or without the use of a yield criterion. A basic assumption 
for this class of constitutive theories is that in the range where inelastic 
strains are present, the total strain rate sjj can be divided into elastic and 
inelastic components which are both nonzero, l.e. 

• ■ t ^ ^ i ^ (1] 

*ij *ij *ij 11 

This equation is applicable for the small strain ease and a similar decompo- 
sition is assumed to hold for the deformation rates in the case of large 
strains. Those are equivalent to strain rates if the strains are small. 

For the small strain case considered here, the elastic strain rate is 
given by the time derivative of Hooke's Law. An important question related to 
Equation (1) is an appropriate definition of *£? . One possibility is to de- 
fine as the total strain rate contribution that is both thermodynamically 
and geometrically irreversible , i.e., non-elastic in all respects. An alter- 
native definition of the incremental plastic strain is the residual strain 
upon loading and unloading from a stress increment. This seems to be the de- 
finition adopted by E. H. Lee in his treatment of large plastic strains, see 
e.g. [14]. Since non-elastic strains are also generated during unloading, 
constitutive equations based on this definition would be different than in the 
former case. Each approach seems possible, but the proper definition and use 
of the term should be indicated and be consistent with the constitutive 

equations. 
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The expression "unified" applied to snch theories is generally taken to 
mean that all aspects of inelastic behavior snch as plastic flow, creep, and 
stress relaxation are included in the function and are particular response 
characteristics for different loading histories. This broad definition of 
"unified" would admit theories with and without a yield criterion and alterna- 
tive specifications of ej?. Separation of the nonrelastic strain rate into 
geometrically reversible (anelastic) and non-reversible components could be a 
convenient procedure and does not detract from the "unified" concept. 

Constitutive theories which are formulated without the use of a yield 
criterion include that of Bodner and his associates [3,4], Walker [2], Miller 
[5], and Crieg, Swearengen and Rhode [6]. Since these models do not contain a 
completely elastic regime, the function that describes the inelastic strain 
rate should have the property that the inelastic strain rate be very small for 
low stress levels. 

For theories with a yield criterion, a^j is identically zero until an 
invariant function of the stress reaches a prescribed value; the function, by 
definition, is independent of strain rate. For stresses at or exceeding the 
yield value. Equation (1) applies and sj? and the stress are function- 
ally related. The fully elastic state, i.e. s^j ■ 0, would apply only for 
stress states less than the rate independent yield value, and loading and un- 
loading paths above that are controlled by the loading conditions through the 
constitutive equations. Theories of this type have been developed by Perzyna 
[15] for the case of isotropic hardening and by Chaboche [7], Robinson [8], 
and Lee and Zaverl [11] for the case of both isotropic and directional harden- 
ing. 


All the unified models are formulated on the basis of internal variables 
which depend on the loading history. The essential features of these unified 
theories are: (1) a flow law which functional form depends on the method of 

treatment of directional (kinematic) hardening, (2) a kinetic equation which 
is the temperature dependent functional relationship between the strain rate 
and stress invariants and includes internal variables, and (3) a set of evolu- 
tion equations for describing the growth of the internal variables. Here, the 
internal variables are used to represent the current resistance to inelastic 
flow of the deformed solid. Two deforming solids with identical values of 
their internal variables would have identical inelastic responses under the 
same imposed stress state. Both the choice and the number of internal vari- 
ables vary with the unified models. • Most of the unified models use two inter- 
nal variables or one variable with two components: one to represent isotropic 
hardening and another to represent directional (kinematic) hardening. In most 
models, the isotropic hardening variable is represented by a scalar quantity, 
either the drag stress (K) or the yield stress (Y) , while directional harden- 
ing is represented by a second order tensor or a scalar function of such a 
tensor. 


Basic Flow Laws 


Four basic forms of the inelastic flow law have been identified. Plas- 
tic incompressibility is always assumed and these flow laws are: 
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where Sjj , ejj , and L ij are the deviatoric, direct and effective stresses, 
respectively. The tensor Qjj represents the "equilibrium stress" which has 
also been referred to as the "bach stress" and the "rest stress." The para- 
meter f is a yield function or a flow potential. It shonld be noted that the 
first three laws can be considered or can be derived from Equation (2d)- if 
they are associated with a flow potential. 

Equation (2a) is the Prandtl-Reuss flow law associated with the von 
Mises yield criterion. However* it can be considered as a basic material 
equation in its own right independently of a yield condition* As such, this 
equation is usually taken to be applicable for proportional loading conditions 
for which isotropic hardening would be appropriate. The equation states that 
the material response (i.e.* the plastic strain rate) to stress is isotropic 
even though Xj could be stress history dependent. Since stress is direc- 
tional. Xj could have a directional character within the context of increment 
tal isotropy and thereby account for induced directional hardening effects. 

Equation (2b) is the flow law obtained by introducing the kinematic 
hardening variable of Prager [Id] into the classical plasticity formulation to 
account for directional hardening (the Bausehinger effect). In this context, 
the term Qj j would represent the new origin of a translating von Mises yield 
surface in deviatoric stress space* and Equation (2b) would be the associated 
flow rule. As before* Equation (2b) can be taken to be a basic material equa- 
tion in a formulation without a yield criterion and the "equilibrium stress" 
tensor Q^j is generally intended to serve the following functions: (a) to 

account for directional hardening (the multi-dimensional Bausehinger effect)* 
and for the non-coaxiality of sj? and Sjj under nonproportional loading his- 
tories (Figure 1); (b) to account for reversed plastic straining effects, e.g. 
reversed creep* relaxation through zero stress* when the effective stress Ejj 
is negative; (c) for theories without a fully elastic range (i.e., a yield 
criterion), to account for low plastic straining within a given range. 

Equation (2c) is the generalized anisotropic form of the Prandtl-Reuss 
flow law which can be rewritten in a 6 D stress and strain rate space to take 
the form. 
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where E a end Tp are related to the usual plastic straia rates sad stresses ia 
a simple manner* see [10] * sad A«0 i*' the 6x6 matrix of coeff icieats. If the 
material is iaitially isotropic sad the law for plastically iadaced direction- 
al hardening does not lead to off diagonal terms* then is initislly and 
remains diagonal. Under these conditions* Equation (2c) is equivalent to 
Equation (2b) since 6 material constants determine the anisotropic flow behav- 
ior. All the flow equations* Equations (2a«b*c)* would be equivalent for the 
case of proportional loading* including cyclic conditions. The real differ- 
ences in those equations would show up for nonproportional loading histories. 


For constitutive theories with a flow potential* both the flow law and 
the growth law of the directional (kinematic) hardening variable Q^j are de- 
rivable from a single flow potential. The associated flow law of a basic form 
of such a flow potential is [8]. 

n-1 

s.? ■ — F 2 (<x , , — Q , , ) for inelastic loading (4a) 

ij 2u 


and ■ 0 


for elastio loading/unloading (4b) 


where F is the von Mises yield function* n and |i are material parameters. The 
conditions for inelastic loading and elastic unloading have been identified 
in [8]. It can be easily seen that (4a) and (2b) are equivalent. In both 
cases* the direction of the inelastic strain rate vector is coaxial with the 
current effective stress vector (see Figure 1). 


Kinetic Equations 

The flow laws* Equations (2a) and (2b) can be squared to give respec- 
tively, 

*1 ’ ® 2 P /J 2 1 1/J <J*> 


H - ® 2 P/J 2 ]1/2 (Si) 

where D 2 P is the second invariant of the plastic strain rate, D 2 P ■ (1/2) SjJ 
8j?* and J 2 end J 2 are the second invariants of the deviatoric stress and 
effective deviatoric stress, respectively, 
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Fundamental to all "unified" viscoplastic formulations based on flow 
laws of the forms listed in Equations (2) is that inelastic deformations are 
governed by a functional relation between D£p and J 2 (or J£) that could in- 
volve load bistory dependent variables. These variables are intended to re- 
present properties of the inelastic state with respect to resistance to plas- 
tic flow. e.g. hardening, and damage. Some funetions that have been suggested 
are the following. 


(a) 

V- 

“o 1 * 


(7a) 
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X ■ 
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and n, m. 

and D 0 1 

ire constants. 

The inelastic strain rate components 

can then 


be obtained as a function of the stress by the use (2a) or (2b) and one of 
Equations (7). Expression (7b) would seem to have some advantage over (7a) 
or (7c) in theories without a yield criterion in that the value of D 2 P is 
almost zero for some range of J 2 regardless of the value of n. In (7b). 

Do is the limiting value of the inelastic strain rate in shear; (7a) and 
(7b) do not contain such a limit. These differences between the kinetic equa- 
tions are illustrated in a normalized plot of log (D 2 P /D Q ) vs X in Figure 2 
for the case of n ■ 3 and m - 1.0. 

In all the preceding equations (7a.b,c) the exponent n influences the 
slope of the D 2 . J 2 relation and therefore has the major influence on strain 
rate sensitivity. That parameter also affects the overall level of stress- 
strain curves although the level also depends on the hardening parameter X. 

Temperature (T) dependence of plastic flow is a first order phenomenon 
comparable to strain rate sensitivity and should appear directly in the ki- 
netic equation. In the case of Equations (7a, b), this can be achieved by 
taking the exponent n to be a function of T, e.g. n ■ ck/T (k is Boltzmann's 
constant and c a material constant) which leads to strong temperature depend- 
ence of the stress parameter X b 3J2/X2 (or 3J2/^) • Numerical results for this 
dependence are shown in Figure 3 for both the power law and exponential kinet- 
ic equations at different non-dimensional ized strain rates. 

The method of including temperature dependence in Equations 7a, b is com- 
parable to an activation energy formulation. Table I lists temperature- 
dependent kinetic equations based on four different functional expressions for 
the aotivation energy. Some of the consequences of the various relations are 
disoussed in [17] . 

Another procedure for including temperature dependence in the kinetic 
equations is to multiply the stress function, the right hand side of equa- 
tions (9) by a temperature function. The temperature factor can again be 
motivated by thermal activation considerations and the Arrhenius expression 
seems to be the reasonable function to use (Table I) . This is the approach 
taken by Miller [5] . 
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Evolutionary Equations fog Internal Variables 


The general framework of the evolutionary equations of internal vari- 
ables is based on the now well-accepted Bailey-Orowan theory [18,19] which 
theorizes inelastic deformation to occur under the actions of two simultan- 
eously competing mechanisms, a hardening process proceeding with deformation 
and a recovery or softening process proceeding with time. The evolution rate 
of an internal variable is then the difference between the hardening rate and 
the recovery rate as given by 

li - h 1 (X i ) - r 1 (X i , T) (8) 

where is the evolution rate of the internal variable X^, and h^ and rj are 
the hardening and the thermal recovery functions, respectively, h^ and rj are 
functions of X±, temperature, T, and the hardening measure, %is either 
or ftp depending on the model. 

(1.) Isotropic Hardening 


The quantity X in Equation (7) is usually interpreted as the iso- 
tropic hardening internal variable and is often referred to as the drag 
stress. Evolutionary equations for the isotropio hardening parameter general- 
ly follow the hardening/ recovery format shown in Equation (8). A comparison 
of these hardening and recovery functions in various unified theories is shown 
in Table II. The rate of isotropic hardening is usually given by a function 
of the hardening variable X, which may saturate to a limiting value, shown as 
Xj in Table II, multiplied by a measure of the hardening rate. Both the in- 
elastic work rate and the effective inelastic strain rate have been proposed 
as the scalar hardening measure. On the other hand, the rate of softening 
or recovery is often taken to be a power function of X and a temperature- 
dependent constant X 0 which value represents the reference state for that par- 
ticular temperature. This recovery model, sometimes credited to Friedel [20], 
theorizes that recovery occurs only when the current internal state exceeds 
the reference state. 

(2.) Directional or Xinematio Hardening 


Probably the main difference in the various unified theories is 
the treatment of directional or kinematic hardening. Differences exist not 
only in the choice of the flow law but also in the evolutionary equations. 
The general framework of these evolutionary equations follows the hardening/ 
recovery formulation represented in Equation (8) with indexes to indicate 
the directions of hardening and recovery. 


h ( V*ij * d( V T) *« • 


W I>v ii 


®(a ,T)i » <9) 


where h 2 » d, and are the hardening, dynamic recovery , and static thermal 
recovery functions, respectively. 0 represents hardening and/or recovery as- 
sociated with the rate of temperature change, ftjj, fly, V^j and are the 
directional indexes of h 2 > d, r2> and 0, respectively. The main differences 
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among the various theories, as summarized in Table III, are in the choices of 
the directional index and the hardening and recovery functions. 

As indicated in Table III, unified models based on the equilibrium 
stress utilize the inelastic strain rate as the directional index for harden- 
ing and contain a "dynamic recovery" term in the hardening function. The pro- 
posed hardening rule is thus similar to the Prager rule [16] in conventional 
plasticity which requires the translation of a yield surface to occur in the 
direction of the plastic strain increment. On the other hand, the evolution- 
ary equation proposed in conjunction with Equation (2a) is based on the direct 
stress as the index for directional hardening [3,4]. This formulation avoids 
the cross-softening effect associated with inelastic strain rate as the index 
and the theory is more compatible with Ziegler's modification [21] of the 
Prager hardening rule. The directional index for "dynamic recovery" is gener- 
ally in the opposite direction of the directional hardening variable Qjj . The 
"dynamic recovery" term is treated in [3,4] as a saturation term in the direc- 
tion of the direct stress but the index has recently been modified to be in 
the direction of - Q^j also [22]. The. unit vector which represents the direc- 
tion cosines of the directional hardening variable is usually taken to be the 
directional index for statio thermal recovery. Recovery always oocurs in the 
opposite direction of the unit vector and tends to reduce the magnitude of the 
directional (kinematic) hardening variable. Most unified theories utilize 
Friedel's recovery model and take zero magnitude of Qjj as the reference 
state. Table III shows that a temperature rate term is also included in the 
theories of Walker and Chaboche. In principle, similar terms could be added 
to the other theories. 

The temperature dependence of the internal variables is also important. 
The experience with the unified models to date indicates that all the material 
constants in the formulations would depend on temperature and must be evalu- 
ated at a number of base temperatures. 

Uniqueness and Stability Criteria 

For stability, unified theories with internal variables must, according 
to Pontor [23], obey the following inequality: 

do £j dej - dXi dii > 0 (10) 

• 

where dojj , dajj, dXi , and dXi represent incremental changes in stress, in- 
elastic strain rate, the current value and the evolution rate of the internal 
variables. The inequality admits classical plastic flow, creep, and stress 
relaxation behavior. It also admits recovery phenomena involving negative 
inelastic work provided that the corresponding changes in the internal vari- 
ables are sufficiently large to make the inequality in Equation (10) remain 
valid* The basic requirement of Equation (10) is that the dissipation rate 
must be nonnegative. 


8 


For a constant internal state* a saall change in <r^j results in a cor- 
responding change in s^ so that [23] 

dffjj dsij >0 , X* - 0 (11) 

The inelastio work inequality is identical to Druoker's postulate [24] in 
classical plasticity that for a stable material flow the plastic work done 
must be nonnegative. For proportional loadings the kinetic equations represen- 
ted in Equation (7a) to Equation (7c) all yield convex "flow potentials" to 
.which the inelastic strain rate vectors are normal. The consequence is that 
the inelastic work is always positive* and unified theories based on Equation 
(7a) to (7c) obey the inelastio work inequality. 

For uniqueness, it appears that the inelastic strain rate must be a 
single-valued function of stress and internal variables. To satisfy the re- 
quirement for stable flow. Equation (10) dictates that stress-strain curves at 
constant strain rate must have positive slope but must decrease with increas- 
ing strain. On the other hand, stress-strain curves at constant plastic 
strain or plastic work must have positive slope, but the slope may either 
increase or decrease with increasing strain rate [23]. 

Most, if not all, of the unified theories listed in Table II satisfy the 
Ponter inequalities and met the uniqueness and stability requirements. The 
stability requirement is, however, not essential for constitutive theory de- 
velopments. Unified theories admit unstable inelastio flow and are generally 
modeled by including softening mechanisms such as thermal softening and con- 
tinuum damage in the evolution and/or kinetic equations. 


NUMERICAL METHODS FOR INTEGRATING 
UNIFIED CON S T I T UTIVE EQUATIONS 


The unified constitutive equations can be characterized as mathematical- 
ly "stiff." That is, in these equations, dependent variables are susceptible 
to large changes from small increments of the independent variables or from 
small time steps. This "stiff" behavior occurs usually with the onset of a 
significant amount of inelastic strains in the loading cycle and is due to the 
generally nonlinear nature of the functional forms that are employed in the 
kinetic equations of these theories. 

A systematic comparison of a variety of approaches for integrating 
unified constitutive equations has been reported by Kumar, Morjaria, and 
Mukherjee [25]. This study concluded that for the constitutive theory of 
Hart, a relatively simple Euler integration method, together with a time 
step control strategy, was optimal when compared with the more sophisticated 
methods. The Valker constitutive theory has been intagrated using the Euler 
single step approach usually without automatic time step control, but rather 
by determining an optimum step size for each problem. Efficiency obtained by 
using this approach has been acceptable and has shown considerable improvement 
over more sophisticated approaches such as higher order Runge— Kutta methods. 
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Tanaka and Millar recently developed a noniterative, self-correcting 
solution (NONSS) method for integrating stiff time-dependent constitutive 
equations [26] . In this approaoh. implicit quantities are removed by Taylor 
expansions of <r, t. and 1 through the incorporation of the integration opera- 
tor a. The method which reduces to the explicit Euler method when a ■ 0 and 
to the implicit Euler method when a > 0 is unconditionally stable for a > 1/2 
and is noniterative. Accuracy is maintained through self-adaptive time con- 
trol and by correcting previous errors st the current step. 

A summary of these various numerical techniques and their applications 
to several unified theories as well as to Norton's law for integrating a uni- 
axial stress-strain curve to a total strain of 1-2% is shown in Table 17. As 
illustrated in Table 17. the explicit Euler is stable when the size of the 
strain increment is kept below 10~4, The size of the strain increment can be 
increased by using an implicit method such as the NGNSS or a-method with a ■ 1 
(implicit Euler). By restricting the comparison to the explicit methods only, 
it appears that there is no substantisl difference between the integrability 
of Walker and Miller theories nor between these unified theories and the 
classical Norton law. The size of the strain increment is. however, somewhat 
sensitive to the values of model constants which describe materiel strain 
rate sensitivity. 


PREDICTIVE AND SIMULATIVE CAPABILITIES OF UNIFIED 
CON S T ITUTIVE THEORIES 


Four of the unified models which have been successfully applied for sim- 
ulating and/or predicting monotonie. cyclic, creep, and stress relaxation be- 
havior are those of Robinson [8]. Wslker [2], Bodner-Partom [3.4]. and Miller 
[5] . Robinson' s model is based on a yield condition and utilizes loading and 
unloading oriteria. while the latter three do not. The kinetic equations com- 
monly used in unified theories without a yield surface or flow potential are 
based on the power-law. exponential, and hyperbolic sine functions; these 
kinetio equations sre represented in Walker, Bodner-Partom and Miller theo- 
ries, respectively. These four theories will be used to illustrate the simu- 
lative and predictive capabilities of the unified theories. 

(1) Monotonic Stress-Strsin Behavior 

All unified theories are capable of reproducing the monotonie stress- 
strain curve. Figure 4 shows an experimental uniaxial tensile stress-strain 
curve of Hastelloy-X deformed at a strain rate of 1.3 x 10~4 sec"l at 922 Z 
and model simulation using Bodner-Partom theory. The computed curve includes 
contributions from both work hardening and thermal recovery. 

(2) Cyclic Stress-Strain (Hysteresis) Behavior 


Bauschinger Effect is represented in most unified theories by a kine- 
matic or directional hardening internal variable. Cyclic hardening, however, 
can be represented by increases in the isotropic hardening variable, the 
directional hardening variable, or both. These different types of cyclic 
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hardening behavior are illustrated in Fignre 5 for Bodner-Partom theory vhich 
does not use an equilibrium stress. The use of different evolutionary equa- 
tions for the equilibrium stress in different regions of stress space allows 
Robinson's model to reproduce rounded hysteresis loops. Examples of cyclic 
saturated hysteresis loops calculated using Robinsdn's model are compared with 
experimental results of 2-l/2Cr-lMo steel in Figure 6 [27] . 

(3) Creep Responses 


Host of the unified models can predict primary and secondary creep re- 
sponses of material subjected to constant load. Steady state creep rates are 
predicted by these unified models to occur when the evolutionary rates of the 
isotropic and/or directional hardening variable vanish as the hardening terms 
are balanced by the thermal recovery. Examples of calculated steady state 
creep rate under constant stress and comparison with experimental data are 
shown in Figure 7 for Bodner-P&rtom' s model. According to the unified theo- 
ries, the steady state creep rate is a function of stress and temperature 
only; it should not depend on the loading histories. This is demonstrated in 
both experimental data and predictions by Miller's model in Figure 8 [28]. 

(4) Stress Relaxation Response 


The behavior of unified constitutive models under stress relaxation is 
analogous to the creep behavior. Under a constant strain condition* the re- 
laxation rate would* again* depend on the ourrent values of the internal vari- 
ables and on the growth laws which describe their changes with time and in- 
elastic deformation. Stress relaxation calculations based on Walker theory is 
compared with experimental data of Hastelloy-X [2] in Figure 9. 

(5) Thetmomechanical Response 


The behavior of unified constitutive theories under thermomechanieal 
cycling depends critically on the change of material constants with tempera- 
ture. In particular, the shape of the predicted thermomechanical loop is 
sensitive to the growth of the kinematic hardening variable (the equilibrium 
stress) with temperature. Walker's model prediction of thermomechanieal loop 
of Hastelloy-X is shown in Figure 10. 

(6) Multiaxial Behavior 


All the unified theories utilize single-valued kinetic equations formu- 
lated in terms of either 3J2/X2 or 3J2/X2. For a constant value of the inter- 
nal variable X and under proportional paths, these kinetic equations predict a 
locus of constant inelastic strain rate invariant in stress space; the shape 
of the predicted "yield surface" or "flow potential" is identical to von Hises 
yield function. For unified models formulated based on the equilibrium 
stress, the size of the "yield surface" is proportional to X, while the center 
of the "yield surface" is at Qjj and translates according to the evolution 
rate of Q^j . 
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Recent studies [22,29], indicate that aaterials exhibit considerably 
sore cyclic hardening when tested under nonproportional paths of combined ten- 
sion and torsion than under proportional paths of tension or torsion only. As 
a result, most if not all, of the constitutive models need to be modified to 
tahe into account the hardening behavior due to out-of-phase loading. 


SUHMAEI AND CONCLUSIONS 


1. A review of more than ten time-temperature dependent elastic- 
viseoplastie constitutive theories indicates that these theories 
differ in the choice of flow law, hinotic equation, and evolu- 
tionary equation of the internal variables. 

2. The unified approach treats all aspects of inelastic deformation 
including plasticity, creep, and stress relaxation using the same 
set of flow law, kinetic equation, and internal variables. 

3. The unified constitutive theories satisfy the uniqueness and 
stability criteria imposed by Drucker's postulate for rate in- 
dependent stable plastic flow and Pouter's inequalities for con- 
stitutive theories based on internal variables. 

4. The unified theories can be formulated either with or without the 
use of a yield criterion. Three basic flow laws are identified in 
theories without a yield criterion. For theories with a yield 
criterion, the associated flow law is derived from the yield 
function or the flow potential. 

5. Three different formulations of the kinetic equations are identified, 
and they include the exponential, power law. and hyperbolic sine 
functions. The exponential formulation gives a limiting inelastic 
strain rate and appears to give better results for high strain rate 
applications. 

5. All three forms of kinetic equations are functions of 3 J 2 /K 2 (or 
3J2/I 2 ) and result in "yield surfaces" and equi-creep rate sur- 
faces which are described by the ^-based von Mises criterion. 

7. The number of internal variables varies among the unified theories. 
Most unified theories use two internal variables, one to repre- 
sent isotropic hardening and one to present kinematic or direc- 
tional hardening. The measure of hardening is either the inelas— 
lastic strain rate or the inelastic work rate. 

8. Directional (kinematic) hardening can be modeled with or without the 
use of an equilibrium stress. The directional index of kinematic 
hardening can be based on either the inelastic strain rate or the 
direct stress. 
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9. Material constants in the unified models are necessarily temperature- 
dependent and reqnired to be evaluated at the temperatures of inter- 
est. There are indications that a temperature rate term is also re- 
quired in the unified theories. 

10. All of the unified theories which are reviewed do not automatically 
predict additional cyclic hardening under nonproportional loading 
paths. Additional texms are needed in the unified theories to i in- 
clude such hardening behavior. 

11. The equilibrinm-stress-based unified theories can describe reverse 
creep and/or reverse stress relaxation behavior without further 
modifications. Unified models which are not based on the equilibrium 
stress would require modification by adding an anelastio term in order 
to take into account these types of behavior. 

12. The unified constitutive equations are stiff but can be integrated using 
either explicit or implicit methods. 
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TABLE I 


FIVE FORMS OF TEMPERATURE- DEPENDENT KINETIC EQUATIONS WITH 
THE CORRESPONDING ACTIVATION ENERGY FUNCTION 


Activation Energy 


Temperature-Dependent Kinetic* Equations 


AH = CAn 





AH = H 0 - Vg(J 2 ) 




H 0 -Vg(J 2 J 

kT 


AH 



D n exp 


H* /, K 2 \ 
kT \3J 2 y 


r j/2 i c/kT 

r / 

,2 \ C/kT l 

* ■ ,T fe] 

Dj? 8 D 0 exp - f 

k) J 

AH 8 Q 

°2 * °0 6x '’[- &i 

f s1nh (~r 


where C, Do* H*, Ho* Q» m, and n are constants; V is the activation 
volume; and k is the Boltzmann's constant. 
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TABLE II 


THE SPECIFIC FORMS OF ISOTROPIC HARDENING AND STATIC THERMAL RECOVERY 
FUNCTIONS USED IN THE SELECTED UNIFIED CONSTITUTIVE THEORIES 

K * h-| (K)M.j - r r (T,K) 


Model 

where M-, = e; e = >/ j £j P e’jJ 

or M-j a W p (Bodner-Partom 1 s Theory) 

Static Thermal 

Hardening Function, b-j(K) Recovery Function, r i(T,K) 

Bodner-Partom 

C-j (K<j - K) 

c 2 (k-k 0 ) p 

Wal ker 

C-j (K-j - K) 

- 

Krieg et al 

C 1 

c 2 (k-k 0 ) p 

Robinson 

C 1 

- 

Chaboche 

CjdC, -K) + fj (i, I, n 1d ) 

- 

Lee and Zaverl 

C-j (K* - K)/ 

- 

Hart 

«i. 


Ghosh 

C,K" <1 

c 2 (k-k 0 ) p 

Miller 

C, CK, -C 4 (s1nh- 1 C 3 |j|)V] 

C 2 [sinh C 3 K 01 ] 11 


where C-j, C 2> Cg, C 4 , Cg, m, p, q, K Q> and K-j are material constants; K* 
Is the saturated value of K; K* is governed by an evolutionary equation 
which is function of e and J 2 . 
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TABLE III 


THE SPECIFIC FUNCTIONS OF ANISOTROPIC HARDENING, DYNAMIC RECOVERY, STATIC THERMAL 
RECOVERY, AND THE TEMPERATURE RATE TERM IN SELECTED UNIFIED CONSTITUTIVE THEORIES 
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TABLE IV 


COMPARISONS OF THE INTEGRATABILITY OF VARIOUS CONSTITUTIVE MODELS 

Comparison Is based on the size of the average strain Increment per step for 
Integrating a uniaxial tensile stress-strain curve to a total strain of 1-2% 


Constitutive 

Model 

Integration Method 

Strain Increment 
Per Step 

Automatic Time or 
Strain Increment 
Control 

Stability 

Reference 

Walker 

Explicit Euler 

10' 5 to 10" 4 

No 

Stable 

Walker [2]} 
Cassentl [30] 


NONSS (a-Method) 






a » ,1 

1 x 10 -4 


Stable 


Bodner-Partom 





Present 


a ■ .5 

1 x 10’ 4 

No 

Stable 

Investigation 


a » 1 

1 x 10- 3 


Stable 



Explicit Euler 

4 x I0* 4 

No 

Uhstable 

Lee et al [31] 

Hiller 

Tayler Series Modified Euler 

4 x 10- 4 

No 

Stable 

Lee et al 


NONSS (a-Method, a-1) 

3 x 10-3* 

Yes 

Stable 

Tanaka [26] 


Explicit Euler 

2 x 10-5* 

Yes 

Stable 

Kumar et al [25] 


Pred 1c tor-Corrector 

3 x 10-5* 

Yes 

Stable 

Kumar et al 

Hart 

High Order Predictor-Corrector 

3 x 10-5* 

Yes 

Stable 

Kumar et al 


Two-Step Adam 

2 x 10-5* 

Yes 

Stable 

Kumar et al 



1 x 10-3 

No 

Unstable 

Lee et al [31 ] 

Norton 

Explicit Euler 

5 x 10- 4 

No 

Stable 

Lee et al 


Taylor Series Modified Euler 

l x 10-3 

No 

Stable 

Lee et al 


* Average strain increment per step = total straln/nuuber of time steps 
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Figure ]. Graphical representation of the basic 
’ flow laws used in the unified consti- 
tutive theories. For theories based 
on an equilibrium stress, the inelas- 
tic strain rate vector is coaxial 
with the effective stress J £fj and nor 
mal to the flow potential f If such a 
concept is used. For theories which 
do not include an equilibrium stress, 
is coaxial with the deviatoric 
stress Sij for both Isotropic and In- 
crementally Isotropic cases but is 
noncoaxial with Sjj for generalized 
anisotropic cases. 


Figure 2. Functional behavior of the kinetic 
equations used In the unified con- 
stitutive theories. The exponential 
formulation in Bodner-Partom's theory 
is seen to give a limiting inelastic 
strain rate of D Q . 
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Figure 3. Functional behavior of temperature- 
dependent kinetic equations utilized 
In Bodner-Partom and Walker theories 
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Figure 4. A mono ton 1c stress-strain curve 
simulated by Bodner-Partom* s 
Model for Hastelloy-X at 1200°F. 

\ 
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Figure 5. Cyclic stress-strain hysteresis loop Figure 6. Stable hysteresis loops of 

simulated by Bodner-Partom's Model 2-l/4Cr-lMo steel for Ac * 

for Hastelloy-X at 1600°F. ±.32% at various strain 

rates at 538°C. The calcu- 
lated curves (solid lines) 
are generated with Robinson's 
Model, from [27]. 
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Figure 7. Steady creep rates as a 
function of stress simu- 
lated by Bodner-Partom's 
Model. 


Figure 8. Miller's Model prediction 
compared with experimental 
data for a creep test with 
a sudden decrease In ap- 
plied stress, from [28]. 
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Figure 9. Negative stress relaxation response of 
Hastelloy-X at 871 °C (1600°F) initiated 
from a steady state hysteresis loop exe- 
cuted at a constant strain rate of 
± 1.35 x 10~3 sec~1 with a strain ampli- 
tude of ± .4%. The calculated curves 
are based on Walker's theory, from [2]. 
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Figure 10. Thermomechanical stress-strain pre- 
diction by Walker's theory. 




